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By high temperature series expansion and exact complete diagonalization the magnetic suscep- 
tibility x{T) and the specific heat C(T) of a two-leg 5* = 1/2 ladder with cyclic (4-spin) exchange 
are computed. Both methods yield convincing results for not too small temperatures. We find that 
a small amount of cyclic exchange influences the thermodynamical properties significantly. Our 
results can serve as reliable basis for an efficient analysis of experimental data 

PACS numbers: 05.10.-a, 75.40.Cx, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 

Recently, it has been pointed out that both for spin 
ladder systems and the parent compounds of high-Tc su- 
perconductors ||l| in addition to the bilinear exchange 
also biquadratic ring exchange terms arc important 
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the modification of the low temperature behavior due to 
this new exchange interaction has been discussed in de- 
tail. However, only little investigations of the impact 
of this new type of interaction on the finite temperature 
properties are available up to now. Therefore, in this 
paper we address the question how the thermodynamic 
properties of two leg spin-1/2 ladders are modified by 
cyclic exchange interactions. In particular, the present 
work provides the appropriate high temperature series 
(HTSE) data and results from exact complete diagonal- 
ization (ED) [|l8j of the magnetic susceptibility and the 
specific heat. Wc expect that such an analysis constitutes 
an important supplement to the study of spin-ladders, 
furnishing additional information about couplings and in- 
teractions. 

The model under study is given by the Hamiltonian 



H — ^J^Si.iS2^i + J|| [Si^iSi^i+i + 82.^82, i+i] 

i 

+ 2 Jcyc [(Si^iSi^i+i)(S2,iS2,i-l-l) (1) 
+ (Sl,iS2,i)(Si^i+lS2,i+l) — (Si^iS2.i+l)(S2^iSi^i_|_l)] 



where Jj_ > and J|| > are the rung and leg couplings, 
respectively; the subscript i denotes the rungs and 1,2 
the two legs. Jcyc > parametrizes the cyclic (4-spin) 
exchange. 

In Sec. |l| we bricfiy sketch the methods wc use. In 
Sec. |l3, results for the magnetic susceptibility x(T) 
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and the specific heat C(T) 
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arc presented for the Hamiltonian in Eq. |^. We also dis- 
cuss representations of our results which can be used to 
pinpoint the parameter regime from experimental data 
quickly. 
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FIG. 1; Two-leg ladder with cyclic (4-spin) exchange 



II. METHODS 

We use ED and the analytic method of HTSE. The 
computer aided calculations for the HTSE yield polyno- 
mials in the coupling parameters with fractions of inte- 
gers as coefficients so that no accuracy is lost. 

Details of our calculations can be found in Ref. |l9|] 
and details of the extrapolation schemes in Refs. plf . 
Further, the data are provided in electronic form so that 
it can be put to use quickly. 

For both quantities xC^) and C(r) the results are pro- 
vided up to order ten in the dimensionless inverse tem- 
perature /3 = J^/T. The first orders for x{T) and C(T) 
are listed in Eqs. [l4a| and 14b, where a;cyc = Jcyc/ J 1. and 
X = J\\/J_i_. Higher orders are available electronically. 

The bare truncated series are not sufficient to describe 
the quantities under study at low values of T. Fade rep- 
resentations are used to enhance decisively the region 
of validity of the HTSE results. Additionally, we include 
low temperature information to enhance the region of 
validity even further. In some cases we can describe the 
susceptibility in the complete temperature regime. This 
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is in particular true if the system is substantially gapped 
because the relevant correlation length remains finite re- 
stricted by the inverse of 'ma'yi{kBT/hv, A/hv), where v 
is a typical velocity of the excitations. 

For a gapped system both C{T) and xi'T) vanish for 
T ^ 0. At finite but small temperatures the deviation 
from zero is exponentially small due to the spin gap. Fur- 
thermore, the leading power in T depends only on the 
dimensionality of the system |Q . The low temperature 
susceptibility for the ladder system can expressed by 



X(T) « T- 



for T < A 



(4a) 



An analogous expression is obtained for C{T) 

C(T)KiT-ie-T' forT<A. (4b) 

Eqs. ^provide additional information to bias the extrap- 
olations. The value for the spin gap is obtained from 
Ref. M. 
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C(r) for various orders of Pade representations and ED re- 
sults from = 12 to = 20. lower plot: Differences between 
HTSE representations. The difference [5,5] - [4,6] cannot be 
discerned in the figure. 



In contrast to our previous extrapolation schemes [pO[ , 
here we employ and extend the method suggested in 
Ref. to extrapolate the HTSE data for the specific 
heat. The main idea is to express the entropy S{T), ob- 
tained from the HTSE data of the specific heat, in the 
new variable e — eo, where eo is the ground state energy 
and e — e(T) is the average energy per site. The tem- 
perature and the specific heat are then derived from the 
entropy as functions of e. The ground state energy is 
obtained from Ref. [Q. To incorporate the low temper- 
ature information from Eq. 1^ the sum rule 
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is considered in the limit T ^ A. Inverting the above 
equation provides an expression T(e — eg). Only an ap- 
proximate solution is possible. Taking the logarithm of 
Eq. 1^ and iterating twice in T leads to 
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with y = e — cq. Combining Eqs. |^ and ^ and the sum 
rule 

S= r ^^^dT' f^T-h-^ for T<A (7) 
Jo ^ 



yields the low temperature behavior of the entropy 
Siy) « -|-ln(\/Aj/)) for y I . 



(8) 



The logarithmic singularity at e = Cq (y = 0) can be 
avoided by extrapolating the function 



G(y) = yd, 
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The value of the gap A is incorporated by requiring 

G{y = 0) = -1/A . (10) 
The entropy is then given by 

Sie) = ie-e.){r$^deJ^'^) (11) 
V Jo e' - eo eo / 

where G is the Fade approximant of G. Wherever possi- 
ble diagonal Fade representations, i.e. the same order in 
numerator and denominator are used. Exceptions to this 
rule will be stated explicitly; they are necessary where 
spurious poles occur. Comparison to ED results show 
that diagonal representations yield the best results, see 
Fig. I The convergence shown in the lower plot is very 
convincing. The plot compares to ED data for system 
sizes up to = 20. It is systematic to the ED calcu- 
lations that for increasing system sizes the results altcr- 
natingly yield an upper (A^ = 12, 16, 20) or a lower 
(A^ = 14, 18) bound on the specific heat for not too low 
temperatures. Therefore, the result for N = oo should be 
in between the results with A^ = 18 and A^ = 20, which is 
convincingly fulfilled by the HTSE extrapolations. The 
shoulders in the low temperature regime at T/ J± < 0.2 
are numerical artifacts. 

It has to be noted that the precise knowledge of the 
ground state energy eo is of particular importance for 
the extrapolation of the HTSE data for the specific heat. 
Even a small uncertainty of half a percent in eo leads 
to significant differences in the specific heat at and be- 
low its maximum. The high temperature part is unaf- 
fected thereof. The ground state energy eo(x,Xcyc) is 
calculated up to order 11 in a; and high order 

series expansion about the limit of isolated rungs. We 
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FIG. 3: Susceptibility (a.) and specific heat (b.) for x = 0.2 
and Xcyc ~ 0.2 for various orders of Dlog-Pade representations 
and ED results for N = 12, 14, 16. HTSE representations 
in b. can not be resolved in the plot. Inset: Phase line with 
A = 0, for details see Ref iiil. 



set r 



^cyc 



/x = const and use standard Dlog-Pade ap- 
proximants on deo{x,rx)/dx which yield highly accurate 
results. 

With the temperature as function of e at hand it is pos- 
sible to represent also the susceptibility as function of e. 
The low temperature behavior from Eq. 4a can also be 
incorporated in the extrapolations. Unfortunately, the 
convergence of the extrapolations investigated is not as 
satisfying as for the specific heat. The low temperature 
regime is underestimated and the extrapolations are very 
sensitive to the order of numerator and denominator in 
the Pade representations. Most diagonal Pade approxi- 
mants are not possible due to spurious poles. Hence we 
refrain from using a xi^) representation. 

The extrapolation of the susceptibility follows Ref . P0[ . 
Basically, the low temperature behavior from Eq. |4a[is 
used to improve the representation. To incorporate the 
low temperature information it is advantageous to map 
the temperature regime to the interval [0, 1] via the sub- 
stitution u ~ + P). All approximants are finite for 
u —I- 1 and the representations are no longer restricted 



to diagonal Dlog-Pade approximants. All representations 
for x{T) are extrapolated with the same order of Dlog- 
Pade approximants to retain a consistent description of 
the HTSE results. We use the [n, 2]-Dlog-Pade repre- 
sentation in u with the only restriction n -\- 2 < 11. 
The representation chosen is checked in the limit of the 
(isotropic) ladder and the Heisenberg chain, where pre- 
cise |6|, ^ or exact results are available. 

To assess the range of validity we investigate various 
orders of Dlog-Pade approximants. They are compared 
to the highest order available. The convergence is very 
satisfying. The orders 8 to 10 differ only by lO^^T/Jx 
from the 11th order for x (not shown). This observation 
does not change significantly for other sets of parameters 
considered. 

Another important check is the comparison to ED re- 
sults. In Fig. data for the susceptibility and the spe- 
cific heat are shown for the gapless point x = 0.2 and 
Xcyc = 0.2, which is exactly known |Q. The inset shows 
the phase line where the gap vanishes. The results are 
obtained by a high order series expansion about the limit 
of isolated rungs. The solid line about the exact point 
shows the highly convergent results. The dotted lines 
give a sketch of the phase line for parameters far away 
from the exact point. The phase diagram is investigated 
in detail in Ref. 1 24 . Left to the phase line the system is 



in a gapped rung singlet phase. On increasing Xcyc the 
gap vanishes linearly in a;cyc and opens again linearly on 
the right of the phase line in a staggered dimer phase 

Mm- ^ 

To derive the low temperature behavior of C and x 
at the gapless point the exactly known triplet dispersion 
uj{q) PC at q — TT is used. A similar analysis as in 
Ref. M yields 



X{T) w ^ for T < 1 



C{T) w Vt for T < 1 



(12a) 
(12b) 



Using the sum rules mentioned above for the specific heat 
leads to an entropy 



S{e) w (e — eo)^^'^ for e — eo small. 



(13) 



The entropy is extrapolated using a Dlog-Pade approxi- 
mant biased to contain the extra information from Eq. |l^. 

The chosen representations of the HTSE results yield 
stable and trustworthy results for the calculated ther- 
modynamical properties. Especially the experimentally 
interesting position and height of the maximum are suf- 
ficiently well described for quantitative predictions. 
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III. RESULTS 

Our aim is to provide results which show the quan- 
titative behavior of the thermodynamic properties and 
in particular the effects of a cyclic spin exchange. With 
the help of computer algebra programs the HTSE results 
can be used easily to determine the model parameters of 
a substance under consideration. Only data of standard 
quantities like the magnetic susceptibility is necessary. 
The occurring ambiguity in determining the model pa- 
rameters ( and see below) from only one quantity like 
x{T) can be resolved by the knowledge of other quanti- 
ties, e.g. the spin gap A or the magnetic specific heat 
C(T) as far as accessible experimentally. 




FIG. 4: x{T) for various values of x for Xcyc ~ 
0,0.02,0.04,0.06,0.08 and 0.1 in ascending order from right 
to left. 




FIG. 5: XmaxTmax versus X for Xcyc = 0, 0.02, 0.04, 0.06, 0.08 
and 0.1 in descending order. The symbols are the calculated 
data points from the HTSE. The filled circle on the right side 
indicates the XmaxTmax value for an isotropic chain which is 
known exactly psl . Horizontal lines show the constant values 
used to rescale x{T) in Fig. ^ 



In Figs. ^ and g we address the information content 
of a measurement of xiT) (cf. [^). Fig. || shows the 
energy-scale independent quantity Xmax?max which is a 
characteristic in experimental measurements. Note that 
for increasing x the differences between the curves for 
various values of Xcyc become smaller, because J|| and 
J± set the changing energy scale, whereas Jcyc stays con- 
stant. In the limit of large x the system approaches two 
independent chains with a decreasing relative interchain 
coupling induced by Jx and Jcyc- 

Once the value of Xmax^max is measured the parame- 
can be read off from Fig. ^ But there 
is still an ambiguity which cannot be resolved by a mea- 
surement of x(T) alone as illustrated in Fig. ^. There the 



A. Susceptibility 

Fig. ^ shows an overview of the magnetic susceptibility 
for various values of the cyclic exchange Xcyc and the leg 
coupling X. The choice of the shown parameter regime 
is taken frorn published values for substances presently 
investigated 0, |l|, |l§ . 

A general behavior for increasing Jcyc at fixed J|| is the 
shift of an increasing Xmax to lower temperatures. This 
effect is induced by the decrease of the whole dispersion, 
i.e. all energies are lowered |, m and the global 1/T 
factor which enhances the value of Xmax- For increasing 
X this effect is weakened. The increasing leg coupling 
provides an additional antiferromagnetic coupling stabi- 
lizing the system against magnetic perturbations. Thus, 
an increasing x stabilizes the system and counteracts an 
increasing Xcyc which destabilizes it. 




FIG. 6: Rescaled susceptibilities for XmaxTmax = 
0.1, 0.095, 0.09, 0.085, 0.08 for the x and Xcyc values shown 
in Fig. I 



rescaled susceptibilities are depicted for the indicated val- 
ues of XmaxJmax in Fig. || (solid horizontal lines). The 
main feature in the susceptibility curves is the maximum. 
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For different sets of parameters for a specific value of 
XmaxTmax the qualitative and quantitative behavior can- 
not be distinguished unless precise measurements in the 
low temperature regime are possible. 

To summarize the latter investigations we state that 
with the present results of ED and HTSE it is difficult 
to ascertain all model parameters from the temperature 
dependence of the susceptibility alone, sec also discus- 
sion in Rcf. [Q. With the knowledge of detailed low 
temperature information this problem can be resolved. 



0.5 




FIG. 7: C{T) for three values for x with Xcyc — 
0, 0.02, 0.04, 0.06, 0.08 and 0.1 in ascending order in direction 
of arrows. Gray lines show [6,4] Fade representations. 



netic specific heat C(T). We are aware, however, that it 
is difficult to extract the magnetic contribution from the 
measured specific heat if the energy scale of the lattice 
vibrations is of the same order as the magnetic couplings. 

Once the magnetic specific heat is known the ambi- 
guity in parameters can be resolved. Fig. |^ shows an 
overview of the magnetic specific heat for x = 0.5, 1, 1.2 
and Xcyc = 0, . . . , 0.1. For increasing leg coupling x and 
fixed Xcyc the position of Cmax shifts to higher temper- 
atures and the height lowers slightly for Xcyc = 0,0.02, 
it stays almost constant for Xcyc = 0.04, 0.06, and it in- 
creases slightly for Xcyc = 0.08,0.1. For increasing cychc 
exchange Xcyc and fixed leg coupling C^ax moves to lower 
temperatures and decreases. This behavior is induced by 
the decreasing overall dispersion, as discussed above. 
IV. CONCLUSIONS 

We investigated the thermodynamic properties of the 
two-leg spin-1/2 ladder with cychc exchange. We used 
two methods, ED and HTSE. The representation of the 
HTSE results were optimized by using Dlog-Pade and 
Fade approximants including the behavior of the consid- 
ered quantities, x{T) and C(T) at low temperatures. 

The results can serve as input for quick and easy data 
analysis to determine the model parameters. Especially 
the experimentally interesting position and height of the 
maxima of x{T) and C(T) are described quantitatively. 

We showed that with the measurement of one quan- 
tity alone, e.g. the magnetic susceptibility it is difficult 
to determine all model parameters. Additional informa- 
tion like the spin gap or the specific heat helps fixing the 
parameters. 



B. Specific heat 

Further information about the magnetic properties of 
substances can be obtained by measuring also the mag- 
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